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ANNUAL REPORT 


I. Project Overview 

Scientific Objectives A 3-year program of research is described that addresses two cases 
of tenuous planetary atmospheres: comets and Io. Funding for a continuation of this project for 
another year has already been tentatively approved. 

Comets In order to understand the physical and chemical processes that produce the 
observed spatial morphology of the cometary coma, it is necessary to analyze observational 
with physically meaningful models. With this in mind, the goals of this project are to complement 
our program of theoretical modeling with observational data analysis regarding the spatial 
distributions of neutral gases and dust in the coma. More specifically the proposed comet 
atmosphere research focusses on three problems. 

• To analyze a set of spatial profiles of CN in comet Halley taken by three observers that 
cover a complete 7.4-day period in April 1986. The object is to specify the time synchronization of 
the parent molecule production at the nucleus (and not just the response of the daughters), and to 
separate the effects of photochemical lifetimes and velocities. 

• To apply our newly -developed dust coma model to various observations. The model is a 
fully 3-D time-dependent treatment of dust particle trajectories which is based on the foundation of 
our H coma model. Realistic physics of light scattering by dust, dust radiation-pressure 
efficiencies, and dust terminal ejection velocities are included. We will address such issues as 
day/night dust production, dust particle size distributions and dust fragmentation-vaporization. 

• To analyze observations of the inner hydrogen coma, which can be optically thick to the 
resonance scattering of Lyman-a radiation, with the newly developed approach that combines a 
spherical radiative transfer model with our Monte Carlo H coma model. Data will include 
published wide field rocket and spacecraft data. 

Ifi The production of atmospheric gases from Io’s volcanoes, the condensation, and 
sublimation of the resulting surface frosts, and the combination of direct (ionization, charge impact 
and charge exchange) and indirect (atmospheric heating and atmospheric sputtering) 
magnetospheric wind-driven processes are important both for shaping the structure of Io's local 
atmosphere and for establishing the conditions for escape of that atmosphere and the ultimate 
injection of heavy ions into Jupiter's plasma torus. In order to establish a theoretical basis for a 
variety of related problems associated with Io's atmosphere and its interaction with the plasma 
torus we propose to apply basic methods and years of experience in treating a variety of 
gasdynamic problems for comets and extended atmosphere studies on Io, and specifically 
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• To study the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and 
plasma given methods and algorithms developed for the study of neutral gas cometary 
atmospheres and the earth's polar wind and plasmasphere. 

II. Progress Report 

The principal parts of the comet portion of the work have gone through the complete cycle of 
planning, work, and publication. One publication in each of the three areas of comet study has 
been published during the second and third years of the project. In these areas the work involved 
application and modification of models which have been developed over more than a decade in 
some cases. For each of these papers, all of the Pi’s contribution is credited to this project alone. 
The Io portion of the work, on the other hand, represents a major new modeling effort. In the 
fourth part of section II. B. we present preliminary results from the Io model, giving the first 2-D 
axisymmetric kinetic atmosphere model. We plan one publication on first results before the end of 
the current project year, and substantial progress toward a major paper on a full model including 
photochemistry, radiational cooling and a complete set of neutral and ion species. 

Cometarv Hydrogen Lvman-a Studies We have applied the coupled radiative transfer/Monte 

Carlo H coma model to the Pioneer Venus Lyman-a observations of comet Halley, as part of a 

collaboration with W.H. Smyth and M.L Marconi (AER, Inc.) in order to verify the approximate 
radiative transfer used in the usual analysis of those wide field scans. As part of this collaboration, 
a paper has been written and published describing the analysis of the spatial profiles contained in 
that data set (Smyth, Marconi, & Combi 1995). The major results are that the hydrogen Monte 
Carlo particle trajectory model (MCPTM), developed several years ago (Combi & Smyth 1988 
a&b), can reproduce the spatial profiles observed by the PVOUVS (Pioneer Venus Orbiter 
Ultraviolet Spectrometer) covering a period from late December, 1985 (r~lAU preperihelion), to 
the time of the spacecraft flybys in early March, 1986 (r~0.85 AU postperihelion). Water 
production rates calculated from the data using the model brings the PVOUVS and IUE water 
production into much better agreement. Furthermore, the results imply that the largest water 
production rates for Halley likely occurred during late January, 1986, one to two weeks before 
perihelion. A reprint of this paper is attached to the Appendix of this report. 

Time-Series Analysis of Cometary Spatial Profiles In the early comprehensive studies of 
spatial profiles of molecular species in comets. Combi & Delsemme (1980b & 1986) warned of the 
influence of the irregular time-variation of the parent gas production rate on the shapes of the 
profiles. Their conclusion was that, provided many spatial profiles taken at different times are 
analyzed, fitted scale lengths (using Haser's model) or lifetimes (using vectorial or Monte Carlo 
models) could be determined in a statistical manner, and the effect of gas production rate could be 
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in a sense be averaged out. This was borne out by the results of Fink, Combi & DiSanti (1991) 
whose Haser model analyses of observations of spatial profiles in comet Halley were generally 
consistent on the averag e with those of other comets, despite the now well-known large amplitude 
periodic variation in the post-perihelion period. 

Combi & Fink (1993) presented a time-dependent model analysis of two sets of pairs of 
spatial profiles of each of C2, CN, NH2 and 0( J D), obtained on March 1 and 2, and April 14 and 
15, 1986. We used the photometry of Schleicher et al. (1990a) as a basis for the time varying 
production rate but allowed for an adjustable phase shift and amplitude enhancement between the 
source at the nucleus and the response of the coma at seen in the photometric aperture. We found 
that spatial profiles, which varied significantly from March 1 to March 2 and from April 14 to April 
15, could be consistently modeled for all four species accounting for time dependence with only 
"average" scale lengths. We were then able to assemble a set spatial profiles of CN constructed 
from observations taken by Dr. Anita Cochran of the University of Texas, and Dr. Rita Schulz of 
the University of Maryland during the same time period which taken together with the first 
observations provide a time-series of CN spatial profiles spanning just over one 7.4-day period 
from April 7 to April 15, 1986. Schulz's observations were particularly interesting because of the 
very large spatial scale covered which allows tracking of the gas production waves (or haloes in 2- 
D) out to 1C)6 km. 

This time sequence of 8-days of CN spatial profiles has been analyzed using the time- 
dependent Monte Carlo model, and can be explained with the same model parameters (scale 
lengths, and phase-lag and amplitude compression of the light curve) as had been used previously 
for only the April 14 and 15 profiles. The most interesting aspect of these results is that the wavy 
features predicted to occur at large distances from the nucleus in the previous paper (that was based 

only on data from within 10^ km from the nucleus), were in fact found in the extended spatial 
profiles from the observations of R. Schulz. A paper describing these results was written, 
submitted and published in the Astrophysical Journal (Combi et al. 1994). A reprint of this paper 
is attached to the Appendix of this report. 

Model Analysis of the Dust Comae of Comets The newly developed dust model was applied 
to study the effect of fragmentation of dust particles on the shapes of isophotes in the dust coma. 
Dust particles when released from the nucleus of a comet are entrained in the expanding gas flow 
created by the vaporization of ices (mainly water ice). Traditional approaches to dusty-gas 
dynamics in the inner comae of comets consider there to be an initial distribution of dust particle 
sizes which do not fragment or evaporate. Dusty-gas dynamics based on the Finson-Probstein 
(1968 a&b) model yields a 1-to-l-to-l correspondence between the size of a dust particle, its 
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terminal velocity owing to gas drag, and its radiation pressure acceleration which creates the 
notable cometary dust tail. 

Using the experience gained in dealing with modeling the H coma (Combi & Smyth 1988b) 
and the newly developed dust coma model, it was clear that the typical elongated shapes of the 

antisun ward isophotes in the dust comae of comets on the scale of > 10 4 km from the nucleus 
requires that the one-to-one-to-one relationship cannot in general be correct. There must be a broad 
range of velocities for any particle size, including small particles having a small velocity but large 
radiation pressure acceleration, in order to explain the elongated shape. A way to create such a 
distribution is if particle fragmentation, or some combination of fragmentation with vaporization, 
routinely occurs within and/or just outside of the dusty-gas acceleration region (i.e. up to several 
hundred km). In this way initially large particles, which are accelerated to fairly slow velocities by 
gas-drag, fragment to form small particles which still move slowly but are subject to a relatively 
large radiation pressure acceleration. Fragmentation has already been suggested as one possible 
interpretation for the flattened gradient in the spatial profiles of dust extracted from Giotto images 
of Comet Halley. Grain vaporization has been suggested as a possible spatially extended source of 
coma gases. The general elongated isophote shapes seen in ground-based images for many years 
represents another possible signature of fragmentation. The results were presented at the 1993 
DPS meeting and published in The Astronomical Journal (Combi 1994). A reprint of this paper is 
attached to the Appendix of this report. 

A Global Kinetic Atmosphere Model of Io The other main topic in tenuous atmosphere 
modeling is to study the global structure of Io's atmosphere its interaction with the Jupiter/Io 
plasma torus, and the large escape rate of neutral gases. The first method we considered for doing 
this was to apply a 2-D axisymmetric adaptive mesh code, which has been acquired from the U. of 
Michigan Aerospace Engineering group in order to construct a lower atmosphere model for Io. 
This is a more advanced and numerically stable model similar to that published by Moreno et al. 
(1991), who had used a fixed polar grid axisymmetric version of the Los Alamos CAVEAT code. 
The Michigan code uses an adaptively-refined mesh of square grids. The second more long-range 
method considered was to make use of the techniques developed for the hybrid gasdynamic-Monte 
Carlo model for comets (Combi 1989; Combi & Smyth 1988 a&b). This would have combined 
the hydrodynamic code for the bottom of the atmosphere with a Monte Carlo particle code for the 
top. In another project we were planning a similar course of action for comet coma modeling. 

Since the conception of this project and the original proposal 3 years ago, a number of 
important and highly relevant studies regarding Io's atmosphere have been published. A number 
of studies now suggest that only a limited area, perhaps 5-20% of Io’s atmosphere is covered by a 

collisionally thick SO 2 atmosphere— i.e. with surface densities perhaps approaching 10^-10^ 
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cm' 3 (Lellouch et al. 1992; Sartoretti, McGrath & Parlesce 1994) in the thickest regions. The 
discovery of the likely molecular ion origin for at least some of Io's fast sodium (Schneider et al. 
1991b) implies an abundance of molecules at high altitudes in the atmosphere. Finally, in work 
performed in a separate project with Dr. William Smyth of AER, Inc., as PI, we (Smyth & Combi 
1995) have constructed a consistent model analysis of sodium "eclipse” observations in the corona 
near Io (surface out to a few Io radii, Schneider et al. 1991a) and the B-cloud emissions (out to 
~100 Io radii). The flux distribution for sodium found in Io's atmosphere is best represented by 
neither a Maxwell-Boltzmann flux distribution nor a classical sputtering distribution. A modified 
sputtering distribution (also called a collisional cascade distribution) provides an excellent and 
consistent picture for producing the entire sodium atmosphere out to 100 Io radii. 

Strobel et al. (1994) have raised the issue that Joule heating may dominate over plasma 
impact heating (Johnson 1989; Popieszalska & Johnson 1992). Both processes are really due to 
the effect of ion-neutral collisions on the atmosphere. In an integrated kinetic model for the 
atmosphere and ionosphere, as being developed here, the point is moot because the actual 
microphysical collision description is explicitly included. 

The significance of this is that Io's atmosphere at large is neither globally collisionally thick 
nor near thermal equilibrium, but that collisions are important Therefore, the application of a 
hydrodynamic model to describe many aspects of Io's global atmosphere is highly suspect, 
although it does retain validity if applied in the right regime, such as in the vicinity of the volcano 
(Moreno et al. 1990). 

There have also been other important developments in work not specifically part of this 
project. In parallel developments, the PI has been working on a fully kinetic Direct Simulation 
Monte Carlo (DSMC) model for cometary atmospheres, and with colleague Dr. Ronald Miller of 
the University of Michigan has been working on various aspects of a plasma hybrid kinetic 
simulation of the Earth's polar wind (Miller & Combi 1994) for the ultimate incoiporation of those 
techniques into an Io model. The DSMC is not the same as the so-called fully kinetic Monte Carlo 
simulations published by Hodges (1990) or more recently by Xie (1994). The latter are iterative 
test particle (TIP) methods, which are only applicable to steady-state problems and which require 
an initial guess for the phase space distribution. A similar ITP model has been applied to the 
hydrogen exosphere (Hodges 1994). On the other hand, a 1-D radial single species hot oxygen 
geocorona exosphere model for the Earth using DSMC techniques has been published recently by 
Shematovich et al. (1994). 

In non-equilibrium gas kinetic theory DSMC models have been used to investigate various 
non-equilibrium aspects of gas flows (such as shock structures, see Bird 1976) for nearly 30 
years. DSMC models are inherently time-dependent, and proceed by following the trajectories of 
many individual molecules, accounting for collisions and other forces. What distinguishes DSMC 
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from ITP is that DSMC is equivalent to a direct finite difference solution of a time-dependent dilute 
gas Boltzmann equation description of the system. The ITP method relys on a more heuristic 
argument of convergence of test particle distribution functions toward the correct distribution 
functions for each species. They have been shown to give reasonable results in some applications 
(Bird 1976). However, ITP is not applicable to non-steady flows, such as would be required for 
understanding system-IU variations in the recent HST observations of S and O near Io, which 
imply variations with the system-EH longitude of Io and when Io passes into and out of eclipse 
(Clarke et al. 1994; Ballester et al. 1994b). The 1-D particle model of Pospieszalska & Johnson 
(1992) employed a method which is intermediate between ITP and DSMC models, but is also 
implicitly steady-state. 

The PI has engaged Dr. Ronald Miller to collaborate on the Io portion of this project for the 
purpose of combining an Earth ionosphere polar wind code with neutral gas algorithms developed 
for the Pi's comet modeling program. Dr. Miller, an Assistant Research Scientist in the Space 
Physics Research Laboratory in collaboration with the PI, is making significant progress toward 
this original long-term scientific goal at a level of effort of only 2 to 3 months per year because of 
his expertise in the development and coding of kinetic models and because the projects from which 
the code is being adapted are already quite mature. Dr. Miller is not new to the Planetary 
Atmospheres program, having worked with Dr. Tamas Gombosi on cometary ionospheres in 
recent years. 

As of this writing the first version of an axisymmetric 2-D code for Io is near completion. A 

1- D spherical version of the kinetic code for a simple hydrostatic equilibrium SC>2 atmosphere has 

been constructed and tested successfully. Figure 1 shows the density and surface wind field for a 

2- D axisymmetric sublimation atmosphere. See the caption for descriptive details. The first 
coupled ion-neutral kinetic model of the plasma-torus/atmosphere interaction is now undergoing 
final testing and debugging and a first paper is expected in the first few months of the new project 
year. 
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Sublimotion Atmosphere — Density (per cubic cm) 



f ig “*: e ** 2D-Axisymmetric Kinetic Model for Io’s Sublimation Atmosphere. Shown are contour plots 
for S0 2 density (above) and horizontal winds (below). The small circle at the bottom is Io’s radius (1820 km) The lower 
horizontal axis is the symmetry axis with the Sun to the left. The radial scale has been stretched so that the large circle 
corresponds to an altitude of 80 km. At low altitudes and near the subsolar point the density and surface winds are similar to 
hydrodynamic models. The kinetic model is valid for high altitudes and subsolar angles >60° where fluid models are invalid. 
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ABSTRACT 

Dust particles when released from the nucleus of a comet are entrained in the expanding gas flow created 
by the vaporization of ices (mainly water ice). Traditional approaches to dusty-gas dynamics in the inner 
comae of comets consider there to be an initial distribution of dust particle sizes which do not fragment or 
evaporate. The standard Finson-Probstein model (and subsequent variations) yields a one-to-one-to-one 
correspondence between the size of a dust particle, its terminal velocity owing to gas drag, and its radiation 
pressure acceleration which creates the notable cometary dust tail. The comparison of a newly developed 
dust coma model shows that the typical elongated shapes of isophotes in the dust comae of comets on the 
scale of >10 1 * * 4 km from the nucleus requires that the one-to-one-to-one relationship between particle size, 
terminal velocity and radiation pressure acceleration cannot in general be correct. There must be a broad 
range of particles including those having a small velocity but large radiation pressure acceleration in order 
to explain the elongated shape. A straightforward way to create such a distribution is if particle 
fragmentation, or some combination of fragmentation with vaporization, routinely occurs within and/or just 
outside of the dusty-gas dynamic acceleration region (i.e., up to several hundred km). In this way initially 
large particles, which are accelerated to fairly slow velocities by gas-drag, fragment to form small particles 
which still move slowly but are subject to a relatively large radiation pressure acceleration. Fragmentation 
has already been suggested as one possible interpretation for the flattened gradient in the spatial profiles of 
dust extracted from Giotto images of Comet Halley. Grain vaporization has been suggested as a possible 
spatially extended source of coma gases. The general elongated isophote shapes seen in ground-based 
images for many years represents another possible signature of fragmentation. 


1. INTRODUCTION 

The standard dust coma treatments originated with the 

pioneering work of Finson & Probstein (1968a,b) which con- 
siders the dust to consist of a population of particles having 
a variation of sizes with a size distribution determined by a 
power law in particle diameter The particles were consid- 
ered to have a constant density, p(a)=po, independent of 
particle size. The light scattering properties which are needed 
not only to describe the observed scattered radiation but also 
to calculate the radiation pressure acceleration on each par- 
ticle size were considered to be simply proportional to the 
geometric cross section of the particle implying that the light 
scattering efficiency, 0^, , and the radiation pressure effi- 
ciency, ? 7 r p, both equal unity and are independent of particle 
size. It is clear from the work of a number of investigators 
that these underlying assumptions are highly oversimplified 
for physical dark absorbing dust particles (e.g., see the ex- 
cellent review by Gnin & Jessberger 1991). 

Jewitt & Luu (1990) have published the results of a 
Monte Carlo dust coma model applied to observations of 
Comet P/Tempel 2. Their model is steady state and spherical. 
Ellis & Neff (1991) have developed a numerical dust model 
for neutral or charged dust for comparison with dust experi- 
ments on Vega 1 and 2 and Giotto. However, their results 
state that the dust particle terminal velocities required to 
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match their model to data are 1.7 times the values implied by 
conventionally dusty-gas dynamics. The resulting dust ve- 
locities would exceed the velocity of the gas that is acceler- 
ating the dust. This seems implausible. 

Baum et al. (1992) have presented a number of models 
compared with CCD images of comets, however their foun- 
tain model upon which they base some of their analysis is 
incorrect owing to a subtle error in their derivation which 
actually had been shown correctly in an old paper by Wallace 
& Miller (1958) as well as in a similar but slightly more 
complicated form by Haser (1966). In their Fig. 5(a), Baum 
et al show a model result for a point source with constant 
radiation pressure which shows isophote contours that, while 
not elongated, appear to be circles drifting tailward. On the 
contrary, the results of Wallace & Miller (1958) shown in 
their Eq. (8) and in the upper left portion their Fig. 3, clearly 
state “that, for any position of the observer whatever, the 
isophotes are circles with center at the nucleus, or arcs of 
circles terminated by the (projected paraboloid) envelope,” 
and “the intensity is inversely proportional to the radius of 
the isophote.” However, the conclusion of Baum et al ., that 
some dust particles seem to decay with scale lengths on the 
order of magnitude of —10 s km, is not dependent on this 
aspect of their model and is consistent with some of the 
Giotto data (see the review by Sekanina 1987). 

The Finson-Probstein method has generally been applied 
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lo wide-held images of I he large-scale dust tail (Finson A 
Probstcin l%<Sh; Sekanina iV Miller 1*1731. The dust coma 
within about l(T km from the nucleus presents a different 
challenge. Fulle (I9S7) introduced a Monte Carlo approach 
to the Finson-Prohstein inversion method which permitted 
relieving some of the restrictive simplifying assumptions. 
Most recently Waniak (1992) has published the results of 
using a similar Monte Carlo approach to the inversion 
method where the calculation was extended to spatially 
anisotropic emission from the nucleus. 

In this paper a dust coma model is presented which is 
based on the Monte Carlo techniques previously applied to 
gas species in cometary atmospheres (Combi & Dclsemme 
1980; Combi & Smyth 1988a) and which is fully time de- 
pendent and three dimensional. The light scattering effi- 
ciency, radiation pressure efficiencies, and dust particle den- 
sities are allowed to take on physically realistic values that 
vary with particle size. Such complexities as anisotropic 
sources, periodically varying production rates, and particle 
fragmentation and decay have been incorporated into the 
model for comparison with data. Anisotropic emission and 
periodically varying production rates have been addressed 
with this model in previous preliminary studies (Combi & 
Fink 1992). The general objective of constructing such a 
model is to compare it with some of the wealth of CCD 
images of the dust coma (as opposed to the dust tail), that are 
available for Comet Halley, in order to explore the possibility 
of seeing visible effects of important physical processes. 

In this paper the shapes of modeled isophote shapes of the 
coma are compared with those detected in the typical field of 
view ( — 10 5 km) afforded by a CCD camera and a fairly 
large telescope. Section 2 of this paper presents the math- 
ematical description of the basic model. Section 3 presents a 
comparison of the standard model with a representative set 
of narrowband continuum images of Comet Halley taken 
during the pre-perihelion period. Section 4 discusses a heu- 
ristic dust fragmentation model which can faithfully repro- 
duce the shapes of the observed isophotes, for different view- 
ing geometries and conditions. And finally. Sec. 5 presents a 
discussion and summary of results. 


2. A GENERAL DUST COMA MODEL 

The basic framework of our dust model is built on the 
basis of the hydrogen coma model discussed in detail in the 
papers by Combi & Smyth ( 1988a, b). It is a fully 3D time- 
dependent Monte Carlo particle trajectory model that builds 
up an entire coma by following the trajectories of many in- 
dividual particles. For the published gas coma model, the 
particles are individual atoms, radicals, and molecules; for 
the dust model the particles are obviously individual dust 
particles. The mathematical treatment for many of the basic 
processes, namely the angular dispersion of initial trajecto- 
ries, time-dependent production rate, heliocentric trajectory 
calculations with solar gravity and radiation pressure, collec- 
tion of particle weights tor the calculation of column densi- 
ties, etc., has already been described in great detail in the 
above papers and will not be discussed further here. What 
will he discussed here are those aspects of the model which 


are unique to the dust model: the physical properties of the 
dust particles, their size distribution and density, terminal 
dusty-gasdynamic velocities, light scattering, and radiation 
pressure efficiency. 

From polarization measurements of Comet Halley (Mukai 
et ai 1987) it is clear that comctary grains have similar op- 
tical properties of a material commonly called astronomical 
silicate: i.e., a dark dielectric. Although that real cometary 
grains are not expected to be spherical but may be quite 
irregular, the adoption of Mie theory which assumes spheri- 
cal particles has proven to be quite useful by other investi- 
gators. The results of Mie calculations by Hoban et ai 
(1989) have been adopted for the light scattering properties 
in our model. 

Hellmich & Schwehm (1984) calculated the radiation 
pressure efficiency as a function of size for dark dielectric 
material, which is believed to be indicative of cometary dust 
particles. The radiation pressure efficiency drops precipi- 
tously for very small particles and peaks for particles whose 
sizes are in the range of the wavelength of the peak of the 
solar radiation. This calculation has been adopted for the 
standard radiation pressure efficiency in the model. 

Probably more uncertain than the question of gross optical 
properties of small dark grains in this problem are their den- 
sities. Although it has been traditional to assume densities of 
1 to 3 gem -3 , it is likely that cometary grains are porous 
aggregates and therefore have a density which decreases with 
increasing radius. Lamy et ai (1987) have suggested a dust 
density of the form 

p(fl) = 2.2- 1 Aa/(a 0 + a ), (1) 


where a — particle radius and a 0 =2 microns. This has been 
chosen as a starting point for the modeling studies presented 
here. They can be easily scaled, although the dependence on 
size for aggregate particles must at least be qualitatively cor- 
rect. Note that the dust denstity function in Eq. (1) has the 
same form as that used by Gombosi (1986) 
[p(a) = 3.0-2.2a/(a 0 + a)] but the dust densities are uni- 
formly somewhat smaller. 

For the differential dust production rate the Hanner-type 
distribution of Gombosi (1986) has been adopted which 
gives 


dz d (a) 4 7r y x k 

da “1“ (1+x 1 ) ’ 


( 2 ) 


where a=dust particle radius, 

z d (a)=i he cumulative dust mass production rate in 
particles per unit time, 

dz d (a)lda= the dust mass production rate for par- 
ticles with radii between a and a + da y 

a {> - a parameter taken to be 2 microns, 
x ~afa {) , 

p tf = density of dust particles of size a, and 
<A, = a normalization constant which is set by the in- 
tegrated dust production rate. 

It is related to the gas production rate with the proportional- 
ity constant such that z d ia m . AX )~Q d - xQ K > and 

z d (a mm ) -- (), where a mw and tf max are the smallest dust par- 
ticle sizes in I he distribution. 
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The value of k is taken to be 1.5 (Gombosi 1986). The 
size distribution function of dust particles produced yields a 
function of the form a (V: 111 where the exponent has a 
slightly smaller magnitude (3.2) for the smallest particles and 
somewhat larger (3.6) for the largest particles, so as to be in 
agreement with various in situ studies (McDonnell et al. 
1991) and with comet tail studies (Waniak 1992, and Fulle 
et al. 1988 for Halley; and Fulle 1987 and 1988 for other 
comets). 

For the particle velocities the best approach is the dusty- 
gas dynamic method initially developed by Probstein (1968) 
and subsequently improved (see Wallis 1982; Marconi & 
Mendis 1983; Gombosi et al. 1985) to include the feedback 
of drag by the dust on the gas flow. Sekanina (Sekanina 
1981; Sekanina & Larson 1983) has developed a simple em- 
pirical method for computing dust particle velocities in a gas 
flow which accounts for the mass loading on the gas flow, the 
gas production rate, and the dust-to-gas mass ratio. Since our 
purpose is to generate an explicitly time-dependent dust 
model, having a fast approximate alternative to having to do 
a quasi-time-dependent treatment by running many indi- 
vidual full dusty-gas dynamic models is of great benefit. In 
any case, any errors intoduced by Sekanina’s approximation 
are, on the whole, much smaller than those introduced by 
the uncertainties in the other physical details and assump- 
tions of the model. 

In Sekanina’s approximation the velocity of dust particles 
with radius, a y is given by 


v a = v t 


0.65 


1 

3 [ 1 + 'l + a(a)- 


VS 0 + 0 . 38 *' 


, 0 . 6 \ 


+ 0.225a(a) l/2 , 


( 3 ) 


where i^ = gas velocity in the dust acceleration region 
(—0.42 km s' 1 ), 

«(a)=[cv g i7 rp (a)]/l«.2y8(a)]. 
c/J=1.83xlO“ 4 gem -2 , 

;^= dust-to-gas mass ratio, 

pMp a a]. 

7 ^= radiation pressure efficiency, 
p fl =d ust particle density [Eq. (1)], 
c/?=0.585X 10” 4 g cm” 2 , 

Z = G^/4tt/? 2 , the vaporization rate per unit surface 
area, 

(2^-gas production rate, and 
/?„ = radius of the nucleus (—6 km for Halley). 
Comparisons of the results of using Eq. (3) with the explicit 
dusty-gas dynamic model results of Gombosi et al. (1986) 
show excellent agreement across the entire particle size dis- 
tribution and for the complete range of gas production rates 
relevant for Halley. 

For the long term secular dependence of the dust produc- 
tion rate, those determined in the photometry of Schleicher 
and co-workers (Schleicher 1993) of r~ 2? for pre-perihelion 
and r 19 for post-perihelion, where r is the heliocentric dis- 
tance of the comet have been adopted. Preliminary results of 
adopting the 7.4 day periodic variation seen in comet Halley 
(Schleicher et al. 1990) and anisotropic coma emission have 


Taiu i I Observational parameters. 


Dale 

am 

Exposure 

(s) 

H 

(AUI 

A 

(AU) 

Phase 

C) 

IVKS Nov V 2529 

300 

1 82 

0.8V 

15.1 

1985 Dec 4.0302 

300 

1 45 

0.66 

36.0 

ms Dec 27.9694 

300 

1.08 

1 06 

54.6 


already been explored using this new dust model (Combi & 
Fink 1992). 


3. COMPARISON OF MODELS WITH DUST COMA IMAGES 

A set of CCD images of Comet P/Halley were made using 
the 61 in. Wyeth telescope at the Oak Ridge Observatory, of 
the Havard-Smithsonian Center for Astrophysics from Sep- 
tember through December 1985. The standard IAU comet 
filter set with wavelengths centered at 3840 A (CN), 4280 A 
(CO*\ 4840 A (continuum), 5160 A (C 2 ), 6840 A (con- 
tinuum) and 7130 A (H 2 0 +) was used and various comet, 
flux-calibrated star, and flatfield exposures were taken. Un- 
fortunately, an unpredictable failure of the camera shutter (a 
sticky shutter) was discovered much later after a consider- 
able effort was made to calibrate the images using standard 
star images and understand the inconsistent results. However, 
the combination of the flatfield and “continuum” images do 
yield images of the dust coma distribution, albeit uncali- 
brated, that are still quite useful. Comparison of several pairs 
of nearly simultaneous 6840 and 4840 A images show that 
there appears to be no major difference in the image gradi- 
ents that could be attributable to the wavelength difference 
(i.e., indicative of some particle size discrimination) so the 
rest of the analysis is based upon the 6840 A images which 
had a much better signal-to-noise ratio because of the better 
sensitivity of the overall system at that wavelength. 

Table I gives the pertinent information for the 6840 A 
images analyzed for this paper. A contour plot display of one 
of the images is shown Fig. 1. Immediately noticeable is that 
the apparent radiation pressure acceleration produces pro- 
gressively more elongated isophotes in the tailward direction 
with increasing distance from the nucleus. This is reminis- 
cent of many other bright comets [see, for example, the iso- 
photometric atlas of Hogner & Richter (1969)]. In this re- 
spect the isophote shapes are similar to the H Lyman-a 
images of a number of comets (Meier et al. 1976; Opal & 
Carruthers 1977; Smyth et al 1991). On the contrary, a com- 
parable model image, shown in Fig. 2, tends to give more 
circular isophotes which terminate rather abruptly in the sun- 
ward direction. The comparison is drawn with the H coma 
images because the reasons for the shape, both the immedi- 
ate phenomenological reason and the ultimate physical cause 
are understood in the case of cometary hydrogen, and thus 
serves as a guide for understanding the similar shape of the 
dust coma. 

Work by Keller and co-workers (Keller & Meier 1976) 
showed that the elongated isophotes seen in the H comae of 
comets required that the H atoms emitted from the inner 
coma possess a broad distribution in speeds away from the 
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Fi(t. 1. A r>N4<l A narrowband CCD image of Comet P Halley on 1985 
December 4 Shown is a contour plot of a 6840 A Halley Watch filtered 
image of Comet P Halley taken with the 60 in. Wveih telescope at the 
Harvard/Smithsoman Oak Ridge Observatory. The image has been corrected 
for flatHcld deviations but is not intensity calibrated. The linear dimension 
of the square image is L.VJHxlO* km and the Sun is clearly to the upper 
left. Note especially the elongated shapes of the contours toward the tail 
direction. 

nucleus. As discussed in Sec. 1 of this paper, Wallace & 
Miller (1958), building on work begun with the original 
fountain model of Eddington (1910), showed that for a 
simple point source emission subject to a constant accelera- 
tion in one direction, the resulting “coma” falls within a 
paraboloid of revolution, and that isophote distribution ap- 
pears to be circles that are centered on the nucleus for all 
viewing geometries . In order to build up a more realistic 
coma, a distribution of particles having a range of radiation 
pressures and/or outflow velocities would be required. Figure 
3 illustrates the appearance of such a single component point 
source model with radiation pressure. In many ways, is it not 
unlike the appearance of the much more complicated dust 
model shown in Fig. 2, despite the addition of a dust particle 
size distribution. Certainly the isophote contours are not 
elongated in the tailward direction as they appear in observed 
images such as Fig. 1. 

In the analogous case of the Lyman-cr coma, the value of 
the radiation pressure acceleration on H atoms is nearly con- 
stant, and only very weakly dependent on the somewhat dif- 
ferent heliocentric velocities of individual atoms. A distribu- 
tion in velocity was found by Keller to be required to 
produce the rather elongated tailward coma seen at large dis- 
tances. This can be thought of as a superposition of many 
Wallace and Miller point source paraboloids. (This discus- 
sion of course ignores the effect of the relative solar orbits 
between various particles and the nucleus which would dis- 
tort the simple point source model, but the basic character is 
still the same. The dust coma model presented here does 



Flo. 2. Standard dust model for Comet Halley on 1985 December 4 A 
surface brightness contour plot is shown for the standard model (no frag 
mentation) as discussed in See. 2. The Sun is up. Note in comparison with 
the data image in Fig. 1 that, although the magnitude of the solar radiation 
pressure induced distortion is approximately correct, the model contours arc 
much more round than the elongated contours in the data in the tailward 
direction. Other data images are similarly elongated. The model consisted of 
1,000.000 dust particles in each of 28 logarithmically distributed size bins 
from 10“ 7 to UF cm. 

include this effect.) A distribution in speeds, formed from the 
sum of three Maxwellian distributions (at 4, 8, and 20 
km s' 1 ), provided a good reproduction of the isophote maps 
of comet Kohoutek (1973 XII). Later, Combi & Smyth 
(1988b) showed that the combination of photodissociation of 
H 2 0 and OH combined with partial collisional thermaliza- 
tion of the H atoms in the inner coma produced a broad 



Fl<i. ?>. Point source mi>del isophote contour map. The appearance of the 
simple point-source mtidel with radiation pressure of Wallace A Miller 
(1958) is illustrated. Note that the isophote contours are circular and cen- 
tered on the position of the nucleus. They end discontinuous^ at the posi- 
tion of the projected paraboloid of revolution whose vertex distance is given 
by r'/2 h. where the value, u . is the outflow speed and the value, h, is the 
radiation pressure acceleration 
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distribution of H atom velocities leaving the inner (< 10* km) 
coma, similar to but more irregular than Keller s Max- 
well ians. This new distribution also provided a good repro- 
duction of the isophote shapes. This time, however, the 
model contained few adjustable parameters, and the speed 
distribution was the straightforward (but complicated) result 
of Ihe actual physical and chemical processes in the coma. 

For dust particles there is a similar problem. Since the 
radiation pressure acceleration and the outflow speed (i.e., 
terminal velocity) of a certain size dust grain depend in a 
similar way on the geometrical cross section and density, one 
finds that there is only a fairly narrow range of the quantity 
v m !2 b, where “y” is the initial radial velocity of a dust par- 
ticle and “ft” is its acceleration, among the particles which 
contribute to most of the total scattered light when integrat- 
ing over the abundance-weighted cross section. If the particle 
velocities are increased or decreased by factors of two, the 
resulting coma stills shows rather circular isophotes, similar 
to the simple Wallace and Miller point source paraboloid 
shape, albeit with a larger or smaller paraboloid envelope. A 
similar effect can be obtained if the size distribution or the 
density is changed. 

Finally, an unusual initial angular distribution of dust 
from the nucleus cannot produce elongated isophotes, except 
possibly for an antisunward “beam 1 ' of dust which just does 
not make physical sense. Various exploratory models have 
been run in connection with this paper and the previous work 
of Combi & Fink (1992), and in none of these cases does one 
obtain tailward-elongated isophotes. Other single velocity 
anisotropic and sunward-throwing point source models ap- 
pear in the paper by Wallace & Miller (1958) and there is 
nothing there to indicate that one could obtain tailward- 
elongated isophotes by some sunward ejection distribution of 
dust. The last argument against a source distribution produc- 
ing the elongation is the fact that it is clearly present in all 
images, with different observing geometry and presumably 
various configurations of active regions on the sunward 
hemisphere. 

In the case of the H coma, radial outflow speeds ranging 
from 1 to more than 20 km s“\ and all having about the 
same radiation pressure acceleration, are required in order to 
produce the elongated isophotes which are observed. In order 
to produce similar elongated isophotes for dust, a similar 
spread in velocities for certain values of the radiation pres- 
sure acceleration would be required. As mentioned earlier, 
provided that dust particles leave the nucleus with certain 
sizes and remain those sizes, dusty-gas dynamics requires 
that for each size (and density) that a single terminal velocity 
will be obtained. One way to account for a variation in speed 
for particles of a certain size is if the acceleration were more 
of a stochastic process (as in molecular collisions). However, 
given the large size of dust particles relative to gas mol- 
ecules, the fluid-drag approach does seem valid. Another 
mechanism which would produce a range in velocities for 
each particle size would be if there were large scale fragmen- 
tation of dust particles within 1000 km from the nucleus. The 
major acceleration occurs within the first several nucleus- 
radii, whereas most ground-based images have pixel resolu- 
tions corresponding to -1000 km. If large dust particles are 
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initially accelerated to rather small velocities by the gas drag 
and then subsequently fragment to smaller sizes, the newly 
created and slowly moving small fragments would be subject 
to relatively large radiation pressure accelerations. This frag- 
mentation process would produce a broad continuum of par- 
ticle velocities for each particle size, and thus produce the 
kind of source needed to account for elongate isophotes. 
Such fragmentation has been suggested by Thomas & Keller 
(1990) and Konno et ai (1993) to explain various aspects of 
the Giotto Halley Multicolor Camera observations of the 
very inner coma of comet Halley. Here, it is suggested that 
there is clear evidence for large scale fragmentation in the 
typical ground-based observatory images of the comae of 
bright comets which have been available for many years. 


4. A HEURST1C DUST FRAGMENTATION MODEL 

To model fragmentation in a rigorous way would require a 
full dusty-gas dynamic treatment, where an original distribu- 
tion of dust particle sizes are initially liberated from the 
nucleus, and the physics of the fragmentation process is pre- 
scribed. Konno et at. (1993) have done something like this in 
their ID steady-state model, where they prescribe a certain 
rate of fragmentation from an initial distribution of particles 
to a set of discrete smaller sizes. However, although thermal 
stress and the dust acceleration process have been suggested 
as likely sources for fragmentation, a real physical descrip- 
tion of the fragmentation process has never been attempted. 

For the purpose of reproducing images, a simple and ad- 
justable parametrized method is required. It is assumed that 
an initial particle of radius, a , is accelerated to the value 
determined by the Sekanina approximation of dusty-gas dy- 
namics. Since the model is already Monte Carlo based, 
Monte Carlo methods have been used to fragment the par- 
ticles. If the particle fragments, then some collection of par- 
ticles such that aj<a will result. If an initial particle of 
radius, a , in the absence of particle vaporization (i.e., pro- 
ducing gas molecules rather fragmentation into smaller dust 
particles), is broken up into equal sized particles of radius, 
fly, then by conservation of mass the number of smaller par- 
ticles will equal 

fl 3 p 

N = (4) 

<>fPf 

If all possible fragment sizes were equally probable then a 
suitable distribution of fragmented particles of radius, ay, 
originating from initial particles of radius, a, could be ob- 
tained using Of=aR i where /?, is a random number on the 
interval from 0 to 1. Since this is not necessarily so, an 
adjustable parameter, p, is introduced such that the new par- 
ticle radius is chosen as 

fly — aRP . (5) 

Such a form can be skewed to favoring smaller or larger 
particles (i.e., more or less fragmentation) by letting p take 
on any value greater than or equal to zero. Zero of course 
would yield the original distribution, that is, no fragmenta- 
tion. 
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Figure 4 shows i he effect of v umlion of l he fragmentation 
parameter. />. on the shapes ol the isophotes lor the same 
model alreaiK shown in Figure 2. The isophotes become 
elongated in the lailward direction just as in the case of the M 
coma. As p gets larger, more smaller particles are produced 
from initially large ones, and the isophotes become more 
elongated in the lailward (anti-sunward) direction. Figure 5 
shows a comparison of three observed images with the 
model having p - 2.0. which provides the best match. Values 
of 1.5 and 2.5 are not as good. Not only are the projected 
isophotes the correct shape but both the general scale of the 
sunward-antisunward asymmetry and the radial profile are 
also reasonably well reproduced. 

The reproduction by the model is certainly not perfect, but 
nor is it yet meant to be. The images of December 4 and 27 
clearlv show effects of sunward-enhanced ejection of dust 
which arc not included in the model. In addition, the outer 
isophote in the lailward direction in the December 27 (to the 
lower left) is more extended in the model than in the data. 
Changing the fragmentation parameter does not effect the 
scale of the sunw nrd-lailward distortion as much as it effects 
the shapes of the isophotes making them more pointed in the 
lailward direction. On the other hand, the inclusion of a 
sunward-enhanced ejection of dust would have the effect of 
suppressing the lailward gradient relative to the sunward gra- 
dient while still preserving the tailward pointed shape indica- 
tive of the effect of fragmentation. Such exploration of the 
very large available parameter space to also investigate spa- 
tial anisotropy, short-term time-dependent variations, as well 
as values of the other dust parameters described in the basic 
model in Sec. 2 will be the focus of future work. 

5. DISCUSSION AND SUMMARY 

Althought the heuristic dust fragmentation distribution is 
certainly not meant to be a unique solution for the detailed 
dust size-velocity distribution, the fact that it can reproduce 
the elongated isophotes, where the standard dust models hav- 
ing the 1-to-l-to-l correspondence between size, velocity, 
and radiation pressure acceleration cannot, does indicate that 
this correspondence cannot in general be correct. That frag- 
mentation occurs in the innermost coma has been suggested 
from the Giotto images (Thomas & Keller 1990 ; Konno 
et ai 1993 ). Thomas and Keller in particular found that each 
original dust particle must fragment into at least 2.7 daughter 
particles. The value implied for the fragmentation parameter, 
p y of 2.0, which best reproduces the shape of the dust iso- 
phote contours, implies a much larger rate of —60 daughter 
particles per parent particle. 

There is, however, a caveat regarding the large rate of 
fragmentation implied here. The optical range CCD images 
show the distribution of particles in the size range as small as 
only about one-tenth of a micron. The many small particles 
detected in Hallev by in situ dust instruments (Sekanina 
19N7) have a verv small radiation scattering efficiency and 
are essentially invisible in optical images. This is the reason 
whv carlv attempts at using the Finson-Probstein method 
yielded a small particle size cutoff in the vicinity of 0 .1-0.3 
/um (see Sekanina ck Miller I *>73). Because of this, the va- 



Fio 4. Variation ol Ihe shape- of the dusl coma with the amount of dusl 
fragmentation. The modeled dust coma for the 19K5 December 4 image ol 
Comet Halley is shown for dust fragmentation parameters of 0 (no fragmen- 
tation). I, 2. T and 4. Models have been run for values at the half-integer 
values. A value ol 20 provides (he best reproduction of isophote contour 
shape. Contour intervals coriesf*md to laelors ol two in surface brightness 
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Nov. 9.2529, 1985 



Dec. 27.9694, 1985 




Fi(i 5. Comparison of the dust fragmentation model with CCD images of Comet P/Hallcy. On the left are shown contour plots of three 6840 A images of 
Comet P/Hallcy on (he dales indicated. On the right arc the corresponding best-fit models with a fragmentation parameter equal, p y equal to 2.0. Note that 
adjacent contour levels arc separated hy a factor of two. Slightly preferential sunward ejection of dust is seen in the observed image on 27 December No 
attempt to fit this asymmetry for the analysis is presented here. Note that in the observations, there is some error in relative positions of the outermost contours 
caused by the uncertainly in the level of the sky background which was taken as the value in one of the sunward most corners of each image where the signal 
was the lowest. Somewhat smaller values improve the agreement in the contour positioning compared with the model. In any case, the shape of the isophote 
contours are not effected by the background uncertainty. Finally, there arc more inner contours in the data images since the data pixel size is 400 X400 whereas 
the model grid resolution is only !00x|00. The linear dimensions of the images from top to bottom are 1 . 840 X JO*, 1.398x10 . and 2.287x|0* km. 
respectively. 


liditv of the fragmentation function cannot be extended to 
much below a tenth of a micron. Whereas the fragmentation 
function is providing a reasonable picture of the velocity/size 
distribution of particles as small as 0.1 gm, it could very 


well be overproducing many much smaller particles, which 
are for all purposes invisible. Unfortunately, images at opti- 
cal wavelengths can provide no constraint on the very small 
particles. Therefore, the fragmentation rate of 60 daughter 
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particles per parcnl particle may be better thought of as an 
upper limit. 

Since this framentation process conserves mass, some- 
thing very close to the original <r v ~~‘ Vft distribution remains, 
in agreement with in situ studies (McDonnell ct al. 1991). 
The initial distribution, before fragmentation, includes very 
large particles (1 m), and so the final distribution is complete 
by any definition below about 0. 1 m, since particles this large 
do not contribute much to the total scattered light in the 
visual dust coma. If an appreciable fraction of the grain mass 
were liberated as a gas in a fragmentation/vaporization pro- 
cess, i.e., if there is grain vaporization in addition to grain 
fragmentation, then it could be treated simply as a loss frac- 
tion. If the loss were independent of initial particle size, then 
somethine close to the same original and final distributions 
would remain. If not, then a very different initial particle size 
distribution might be required, in order to yield the “ob- 
served” distribution as a by-product. 

The CCD images of Baum ct al. (1992) provided evi- 
dence for grain fading, i.e., a departure from 1/p in the iso- 
photes. in 10 out of 14 comets they observed. This fading 
was explained either by a decrease in the albedo or the size 
(more likely) of particles with time. Combi & Fink (1992) 
reported that, at least for the case of P/Halley, any departures 
from 1/p were understood in terms of the well-known time 
varying production rate and did not require grain fading. 
Fading of the very small, suboptical grains could of course 
occur without effecting the appearance of optical images. 

The fragmentation of dust particles (Thomas & Keller 
1990; Konno et al 1993) has the opposite effect of fading on 
the radial profile, since as particles break apart, the total sur- 
face area of all the small particles is larger than that of the 
original large particles. On the contrary, fading or vaporiza- 
tion of large grains into small grains via the production of 
gas produces a smaller total surface area in the resultant 
small particles since the number of particles is conserved. 
Although, the analysis of the ground-based images in this 
paper cannot per se distinguish between fragmentation and 
vaporization as the primary source of the required slow and 
small grains, the appearance of the flatter than 1/p distribu- 
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lion in the Giotto HMC images in the immediate vicinity of 
the nucleus would seem to favor fragmentation being domi- 
nant over vaporization. Since dust-to-gas mass ratios are 
thought to be of order unity, certainly 10% of the dust mass 
could be losl to vaporization, and explain an additional ex- 
tended source of gas, whereas 90% of initially large dust 
grains could fragment into many smaller dust grains, and still 
substantially increase the total scattering surface are, thereby 
producing the required flat radial profile. 

Fragmentation and vaporization of dust particles in the 
region immediately suirounding the nucleus have been sug- 
gested to be a possible mechanism for expanding the sub- 
sonic flow region, isotropicizing the day-side source, serving 
as an extended source of gas, and being responsible for the 
surface breeze that sends dust directly to the night-side of the 
coma very near the nucleus (Keller et al 1990). Clearly the 
combination of fragmentation and vaporization of dust in the 
innermost comae of comets may be one of the most impor- 
tant mechanisms for explaining a number of poorly under- 
stood phenomena in comets. 

Further lines of important study in this area are the more 
careful analysis of more dust coma and tail images of comets 
with a fragmentation model, and the two- and three- 
dimensional modeling of dusty-gas flows with dust fragmen- 
tation and vaporization. Future applications of the dust 
model presented here will include dust coma and tail image 
analyses with anisotropic ejection models [done in prelimi- 
nary fashion already by Combi & Fink (1992)], as well as 
more careful study of the constraints (or lack thereof) pre- 
sented by the small (sub-optical) grains. 

The unpublished narrowband CCD images of Comet 
P/Halley shown here were taken and reduced by Dr. Richard 
E. McCrosky (Harvard/Smithsonian), who in turn thanks 
GTE Research Laboratories for instrumentation support to 
Oak Ridge Observatory. The critical reading by an anony- 
mous referee and helpful suggestions is gratefully acknowl- 
edged. Support for this work was provided by Grant No. 
NASW-1908 and NAGW-3417 from the Planetary Atmo- 
spheres program. 
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ABSTRACT 

CN profiles in comet P/Halley were constructed from observations taken at three observatories during an 8 
day period in April 1986. These data provide a time series of CN spatial profiles spanning just over one 7.37 
day period from 1986 April 7 to April 15 and sample distances from the nucleus from just over 10 3 km to 10 6 
km. The effect of the 7.37 day periodic variation on the CN distribution in P/Halley has been examined by 
using the time-dependent model applied earlier to a subset of the data. Because of the large spatial scale of the 
data on April 7, 8, and 9 (- 10 6 km), and the corresponding transport time in the coma, information present 
in the spatial profiles regarding the gas production rate actually covers nearly two full periods. These spatially 
extended profiles clearly show the wavy structures outside 10 5 km. Such structures were predicted in a pre- 
vious analysis (Combi & Fink 1993) that was based solely on the photometric light curve and on profiles 
which only extended to distances less than 10 5 km. We are now able to reproduce the highly variable Halley 
CN profiles using the same model as before, with the same standard scale lengths, phase lag, and amplitude 
correction for the variation in gas production rate. 

Subject heading: comets: individual (Halley) 

1. INTRODUCTION 

In the earlier comprehensive studies of spatial profiles of the 
molecular species, CN and C 2 , in comets, Combi & Delsemme 
(1980, 1986) warned of the influence of the irregular time varia- 
tion of the parent gas production rate on the shapes of the 
profiles and its influence on dissociation scale lengths extracted 
from those profiles. Their conclusion was that, provided many 
spatial profiles taken at different times are included and that 
fitted scale lengths (using Haser’s model) or lifetimes (using 
vectorial or Monte Carlo models) are treated in a statistical 
manner, the effect of gas producton rate could be, in a sense, 
averaged out. This was homed out by the results of Fink, 

Combi, & DiSanti (1991), whose Haser model analyses of 
observations of spatial profiles in comet Halley were generally 
consistent on the average with those of other comets, despite 
the now well-known large amplitude periodic variation in the 
postperihelion period (Schleicher et al. 1990). 

Although the detailed nature of the compound rotational 
motion of the nucleus of Halley’s comet is not known for 
certain (Samarasinha & A’Hearn 1991 ; Belton et al. 1992), it is 
clear that throughout its active period, and especially during 
the postperihelion period, the rotational motion caused a 
variation in gas and dust production with a period between 7.3 
and 7.6 days (Millis & Schleicher 1986; Schleicher et al. 1990). 

The effect of the periodically varying gas production rate was 
demonstrated by Combi & Fink (1993) in spatial profiles of 
CN, C 2 , NH 2 and 0( 1 D) taken on 1986 March 1 and 2 and 
April 14 and 15 by a model based on the photometric light 
curve of Schleicher et al. (1990). All profiles covered radial 


distances of less than 10 s km, and so each only sampled just 
over 1 day of cometary activity out of the total 7.37 day period. 
Furthermore, the profiles predicted a continuation of wavy 
structures in the profiles which should continue out to large 
distances (10 6 km) and have detectable changes with time. 

In this paper we present an analysis of spatial profiles of CN 
in comet P/Halley constructed from observations that were 
taken at three observatories and covering 8 days, from 1986 
April 7 through April 15. Because of the large spatial scale 
(~ 10 6 km) of the data from the early part of this period, and 
given the transport times from the nucleus to these distances, 
the data actually contain information about two full 7,37 day 
periods. The model of Combi & Fink (1993) with the same 
parameters for the scale lengths and the phase lag and ampli- 
tude correction to the photometric light curve can reproduce 
the structure in all the profiles. This paper first describes that 
observational data, then the model analysis and results, fo 1 
lowed lastly by a discussion and conclusions. 

Z OBSERVATIONS 

Spatial profiles of the CN (0-0) band of the violet system 
and the (1-0) band of the red system were collected from three 
separate observing campaigns and are summarized along with 
the Comet Halley’s geometry in Table I. Data were obtained 
from ESO/LaSilla, McDonald Observatory of the University 
of Texas, and the Lunar and Planetary Laboratory of the Uni- 
versity of Arizona. For simplicity, the three data sets will be 
referred to as ESO, MCD, and LPL, respectively, in the 
remainder of this paper. 
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TABLE I 

Spatial Profiles of CN in Comet P/Halley 


Dale R A Profile 

(UT 1986) (AU) (AU) Site* Orientation 


April 7.09* 1.27 0.44 ESO Perpendicular 

April 8.09 1.29 0.43 ESO Perpendicular 

April 9.14 1.30 0.42 ESO Perpendicular 

April 10.40 1.32 0.42 MCD Sun/tail 

April 14.3 1.38 0.43 LPL Sun/tail 

April 15.3 1.39 0.44 LPL Sun/tail 


Notes. — R is the heliocentric distance in AU, A is the geocentric 
distance in AU, Sun/tail refers to the average of spatial profile along the 
sunward and tailward directions of the comet-Sun line, and perpendicu- 
lar refers to the spatial profile perpendicular to projection of comet-Sun 
line. 

* ESO is ESO/LaSiIla, MCD is McDonald Observatory, and LPL is 
University of Arizona. 

The data for April 7, 8, and 9 (ESO) were constructed from 
wide-fieid filtered photographic images of the comet taken by 
the Bochum Hailey Monitoring Program at ESO. The data 
were taken with an f/4 Lichtenknecker Flat-Field Camera with 
a focal length of 760 mm and a corresponding field of view of 
1?8 by 2?7. They were digitized with a PDS microdensitometer 
and calibrated for relative intensity using spot sensitometry. 
Other details regarding these data and their reduction have 
already been published (Celnik et al. 1988; Schulz & Schlosser 
1990; Schulz 1992). 

The profile from April 10 was constructed from long-slit 
spatial-spectral images obtained with a 512 x 320 pixel RCA 
CCD with the Large Cassegrain Spectrograph on the 2.7 m 
telescope MCD. The original spectra covered wavelengths 
from 3650-6800 A. The slit was 2?1 wide for most of the obser- 
vations. A pixel corresponded to 173, but observations in the 
outer coma were binned to 2*6 to improve the signal-to-noise 
ratio. Seven separate observations were made, stepping the slit 
through the coma by computer control of the telescope. The 
general data set and more details regarding the observations 
have already been published elsewhere (Cochran & Cochran 
1990). 

We include here for completeness the two CN profiles of 
April 14 and 15 (LPLX already analyzed in the previous paper 
by Combi & Fink (1993). These data were obtained from long- 
slit CCD spectra covering wavelengths of 5000-10,000 A that 
were taken with the 154 cm Catalina Site Telescope of the 
University of Arizona Observatories. Contrary to the above 
data, these correspond to the (1-0) band of the CN red system 
at 9200 A. The observations, data reduction, and construction 
of spatial profiles have already been presented in detail by Fink 
et al. (1991). Their paper points out the agreement between 
spatial profiles taken in mid-December of the (1-0) red system 
with those of the (0-0) violet system which appear in the paper 
of Cochran & Barker (1986). Therefore, it is clear that the red 
and violet systems can both be used to construct spatial pro- 
files, although the violet system does provide more signal and 
far less background from dust-continuum and night-sky con- 
tamination. 

Although an absolute calibration for the spatially integrated 
long-slit spectra is now available (Fink 1994), none was used 
for the spatial profiles (Fink et al. 1991). The MCD and ESO 
data had their own separate calibrations. However, inter- 
comparison between the data revealed that some of the absol- 
ute calibrations were suspect, especially those obtained at very 


large zenith angles (large air mass). Because of this, we 
obtained CN column densities from Schleicher (1992) which 
are updates and more complete than those originallv Dublished 
by Millis & Schleicher ( 1 986). ^ 

A confusion zone close to the nucleus is present in the spa- 
tially resolved data and is caused by the unpredictability of 
guiding and seeing with slit or aperture size and averaged over 
the source region of the coma where brightness gradients are 
large. Because of this, we chose to calibrate the spatial profiles 
by using the difference between pairs of small and large aper- 
ture photometric observations. The two aperture sizes were 
typically about 80' and 150* in diameter. The difference in the 
column abundances between two such apertures yields the 
total abundance in an annulus with inner radius equal to the 
small aperture and outer radius equal to the larger one. Each 
radial profile was then treated as if it were in arbitrary units 
and integrated numerically around an annulus of the same size, 
assuming the distribution was circular. The ratio between the 
photometric abundances and the spatial profile integral 
yielded a simple calibration factor, 5, for each spatial profile, 
which is given as 

5 _ m Al -m As 

2k pB(p)dp ’ 

where M At is the number of molecules in larger aperture, 
is the number of molecules in smaller aperture, p is the project- 
ed distance from the nucleus, B(p) is the brightness as a func- 
tion of p, and M Al and M As are obtained from the photometry 
(Schleicher 1992). Since the IHW photometry net collected 
data from around the world, there were high-quality photo- 
metric observations obtained within minutes of each of the 
spatial profiles. The photometry results used in our absolute 
calibration is listed in Table 1 

When this scaling was done for the already calibrated MCD 
and ESO data, values ranged from 1 to 2 even after accounting 


TABLE 2 

Results of CN Photometry Used in Absolute Calibration 
of Spatial Profiles 


Data 

(UT 1986) 

\otf 

(km) 

k*/ 

(«*»*■') 

log AT 

April 7.092 

4.09 

-12.473 

30.62 

April 7.098 

4.39 


31.13 

April 8.079 

4.08 

- 12.472 

30.72 

April 8.112 

4.38 


31.16 

April 8.283 

4.08 

-11473 

30.75 

April 8.290 

4.38 


31.20 

April 9.117 

4.07 

- IZ473 

30.77 

April 9140 

4.37 


31.23 

April 10.395 

4.07 

— 1 1474 

30.51 

April 10.401 

4.37 


31.01 

April 14.352 

4.08 

-12.475 

30.54 

April 14.374 

4.38 


31.05 

April 15.296 

4.09 

-1Z476 

30.60 

April 15.305 

4.39 


31.01 


Note. Informalion in this tabic was provided by Schleicher 
1993 . 

• The term /> is the radius of the photometric aperture. 

b The ^-factor is the radiative emission rate for CN(0-0) 
reduced to I AU in ergs s " 1 . 

* M is the number of CN radicals within the aperture. 
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for differences in adopted (/-factors. In theory, all values should 
have been unity. However, the use of standard observatory 
extinction curves at high air masses, rather than measured 
ones, and observing at unusually high air mass made the inter- 
nal calibrations very uncertain in two cases. The consistent 
external calibration of all spatial profiles using the annuli from 
the photometric observations yielded excellent results. It 
should be mentioned that the analysis presented here is sensi- 
tive to the relative calibration between the spatial profile data 
rather than the absolute calibration per se. It was clear that 
there was a relative calibration problem between the indepen- 
dently calibrated spatial profile data sets which were alleviated 
by the use of the photometric data set. 

3. MODEL ANALYSIS AND RESULTS 

The results of Combi & Fink (1993) showed that a time- 
dependent model constructed with Monte Carlo methods and 
the photometric light curves of Schleicher et al. (1990) could 
reproduce pairs of profiles for four molecular species in comet 
Halley on March 1 and 2 and April 14 and 15. The standard 
March and April light curves of Schleicher et al. (1990) had to 
be shifted in time, implying an expected phase lag between the 
activity at the comet nucleus and observed changes in the large 
photometric aperture. The light curve also had to be 
“ decompressed ” because the amplitude seen in the photo- 
metry is compressed both by aperture integration and by the 
extended source nature of producing daughter species, as com- 
pared with the amplitude of the variation of the parent gas 
production source. For April data, a phase lag of 6 hr and an 
amplitude correction of 27% was found to be required for all 
four species in question, CN, C 2 , NH 2 and O^D). 

For the analysis in this paper we simply tried the same 
model parameters for the times (i.e., the phases of the periodic 
variation) that each of the new observations were made. In the 
previous analysis of 2 adjacent days of data, the change in 
production rate owing to the long-term variation with helio- 
centric distance (r -2 postperihelion), that Schleicher et al. 
(1990) had divided out of their results before finding the 
periodic variation, yielded only a small shift ( — 2%) between 
data from 2 adjacent days. For the 8 days of data presented 
here, this long-term variation yields a 25% ramp down in the 
production rate from April 7 to April 15. This had to be taken 
into account for this analysis, the results of which are shown in 
Figure 1. 

The points are the data, and the lines are the model results. 
The absolute scale of the column density corresponds to the 
top profile from April 7. All of the others, both model and data, 
are shifted down by factors of 3 in column density for clear 
examination. The results for April 14 and 15 are the same as 
those described earlier by Combi & Fink (1993), except for the 
small scaling factor which accounts for the long-term pro- 
duction rate variation. Clearly, the wavy features seen at large 
distances, in the range from 10 5 to 10 6 km from the nucleus, as 
predicted in the published models, are evident in the data from 
April 7, 8, and 9. It should be stressed that there was no inde- 
pendent adjustment for each profile. The relative positioning of 
all model profiles with all observed profiles was locked as a 
group. 

4. DISCUSSION 

The general shapes and variations of the observed profiles 
arc well reproduced by the model, but are not perfect. There 
arc a number of reasons why some deviations should be 



10 3 10 4 10 s 10 6 

Radius (km) 

Fig. Comparison of 8 days of CN spatial profiles with the lime- 
dependent model. The points give the observed radial profiles covering the 
period from April 7 through April 15 [from top to bottom ). The lines corre- 
spond to the model of Combi A Fink (1993). The column density scale corre- 
sponds to the first (April 7.09) profile and each subsequent one (both data and 
model) is shifted downward by a factor of 3 below the one above so they can be 
distinguished. 

expected. The model is a Haser-like time-dependent model, so 
all CN radicals have radially outflowing velocities. In reality 
there will be some dispersion caused by the 44 vectorial ” ejec- 
tion of radicals when produced from the parent molecules (or 
grains). At times closer to perihelion, when coma densities were 
larger and dissociation lifetimes shorter, more radicals would 
have been produced well within the collision zone where much 
of this excess velocity from dissociation would have been col- 
lisionally quenched (see, for example, the OH models results of 
Combi, Bos, & Smyth 1993). The fact that radicals will then 
have some range of velocities as they move away from the 
nucleus then implies that the peaks and troughs of the varia- 
tion (as in a wave propagation) should be compressed at larger 
distances from the nucleus. The model does not include this 
effect, so this behavior is apparent in the model-data compari- 
son. 

The CN radicals seen at distances of 10 6 km from the 
nucleus on April 7 were produced as part of a parent molecules 
and/or grains 12 days earlier, near March 26 when the helio- 
centric distance was 1.1 AU, as compared with 1.3 AU in mid- 
April. Furthermore, the production rates and corresponding 
coma outflow speeds were higher at earlier times. There is 
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evidence in the model-data comparison I hat (he wave struc- 
tures have propagated a bit farther out in the data compared 
with the model. This could be explained by this expected varia- 
tion in outflow speed. It is also clear from the original photo- 
metric data base that every 7.4 day period of activity was not 
identical to every other (Schleicher et at. 1990), so some real 
differences are probably present since the model is based on the 
shape of an average light curve. 

Finally, it should be mentioned that the model is obviously 
spherically symmetric and does not make an attempt to dis- 
tinguish between CN jets (AHearn et al. 1986) and CN back- 
ground (Klavetter & A’Hearn 1992, 1994). Deviations from 
circularity were apparent in some of the spatial profiles in 
opposite directions, but these have been averaged out by 
averaging radial profiles on opposite sides of the nucleus. 

In order to analyze wide-field spatial profiles or images 
covering distances out to 10 6 km, it will be necessary to make 
use of a completely lime-dependent calculation which accounts 
for both long- and short-term variations in production rate 
and for long-term variations in photochemical lifetimes and 
radiation pressure acceleration. This would be akin to a heavy 
radical version of the hydrogen coma model constructed to 
analyze wide-field Lya images of comets (Combi & Smyth 
1988). Although the typical spatial scales of those images are 
an order of magnitude larger, the H velocities are also an order 
of magnitude larger, so the problems are quite scalable. 

As mentioned above, the one-dimensional profile data pre- 
sented here do not contain any information regarding the 
angular distribution of material ejected from localized areas on 
the rotating nucleus. This particular isotropic model analysis 
of these data is not sensitive to any changes in the spatial 
orientation of active areas and therefore cannot provide any 
discrimination between alternative nucleus models. It should 
be noted, however, that the 7 day periodically varying light 
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curve (Schleicher et al. 1990), which is used as the basis of our 
time-dependent model, is also one of the important data sets 
contributing to the phasing and period of the nucleus models 
(Samarasinha & A'Hearn 1991 ; Belton ct al. 1992). 

The results presented here do indicate that it may be pos- 
sible to use a more complicated model to analyze a large set of 
two-dimensional images in order to test the complex rotational 
motions in the alternative models of the nucleus. This would 
not be without some further model assumptions, though, since 
the flow physics responsible for producing the coma that is 
somewhat isotropicized from what appear to be obviously 
localized sources on the sunward side of the nucleus (Keller, 
Marconi, & Thomas 1990) is not yet understood. Fragmenta- 
tion and continued vaporization of dust have been suggested 
as a possible source for day-to-night side breezes which could 
inject some dust directly into the night-side hemisphere and 
would also be an obvious explanation for the somewhat iso- 
tropicized (Klavetter & A'Hearn 1994) gas flow in the coma. 

Finally, since the use of simple time-dependent models with 
standard average parameters (scale lengths and lifetimes, 
velocities) can reproduce spatial profiles which are otherwise 
shaped by variable production rates, as in the periodic cases of 
Halley, P/Swift-Tuttle, or Levy (1990c), then a single or limited 
set of spatial profiles could be used to generate the light curve 
of cometary activity with the appropriate inversion technique. 
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The Pioneer Venus Orbiter Ultraviolet Spectrometer measure- 
ments of the Lyman -or intensity of atomic hydrogen excited by 
solar resonance scattering in the coma of Comet P/Halley ac- 
quired from December 28, 1985, to January 6, 1986, and from 
January 31, 1986, to March 6, 1986, are simulated with the Monte 
Carlo Particle Trajectory Model corrected for optical depth ef- 
fects. Spatially detailed comparisons between data and model 
show excellent agreement and are used to infer that the highest 
cometary activity may not be at perihelion, but about 2 Vi weeks 
before. An improved set of H 2 0 production rates is presented for 
the period of time that the spectrometer was observing and is 
found to be consistent with the rates from other types of measure- 
ments. The apparent discrepancy between Stewart (1987, Astron, 
Astrophys. 187, 369-374) in early March and International Ultra- 
violet Explorer OH derived rates is resolved. The problem with the 
conversion of 18-cm OH radio brightnesses to H 2 0 production 
rates is also discussed, e 1995 Ac*kmk p b*. 


1. INTRODUCTION 

The 1986 passage of Comet Halley through perihelion 

was received with an unprecedented observational effort. 

Between September 1985 and June 1986, Comet Halley 

was the target of long-term systematic observations by 

multiple groundbased and orbit based facilities as well as 
in situ measurements by five flyby missions. Observa- 
tions were carried out from the radio to the ultraviolet. 
Part of the vigorous program was the Pioneer Venus Or- 
biter Ultraviolet Spectrometer (PVOUVS) observations 
of solar Lyman-a resonantly scattered from the hydrogen 
in the coma of Halley. Due to the proximity of Halley to 
the Sun in the terrestrial sky, only PVOUVS Lyman-a 
and OH 18-cm radio measurements of Halley could be 
taken near perihelion. Since the OH brightness radio 
measurements for large production rates such as oc- 
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curred for Halley near perihelion are not easily inter- 
preted, the PVOUVS data attain particular importance. 
As a result, although Stewart (1987) already analyzed the 
PVOUVS data, here the data are reexamined with a sig- 
nificantly improved approach. 

Stewart (1987) studied the data using a Haser model, 
neglecting the optical thickness of the hydrogen coma in 
Lyman-a and utilizing a brightness which represented a 
spatial average centered on the brightest point of the 
data. However, as shown in Smyth et al, (1991) and also 
below in this paper, optical depth corrections are quite 
substantial near the brightest point. In addition, the spa- 
tial averaging of the data by Stewart (1987) eliminated 
valuable spatial information, and the scale of averaging 
by early March 1986 was so large (5 x 10 1 * * * * 6 km by 8 x 10 6 
km) that temporal resolution was also compromised. 
Nevertheless, Stewart’s H 2 0 production rates seem rea- 
sonable overall, with the only apparent discrepancy oc- 
curring in early March when the International Ultraviolet 
Explorer (IUE) OH measurements (Combi et al . 1993) 
indicate a much lower rate. 

In this article the Monte Carlo Particle Trajectory 
Model (MCPTM), which has been used very successfully 
to study various other data sets such as IUE measure- 
ments for OH (Combi et al . 1993) and H (Combi and 
Feldman 1992, 1993) and groundbased measurements for 
Ha (Smyth et al, 1993) as well as 0('D) and NH 2 (Smyth 
et al. 1995), is adopted. As discussed below, the MCPTM 
attempts to treat properly the dynamics and chemistry of 
cometary hydrogen. The MCPTM is additionally cor- 
rected for optical thickness effects, and as shown in 
Smyth et al, (1991) as well as here, it is able to accurately 
simulate the spatial brightness distribution from 
PVOUVS. Consequently, it is possible to extract infor- 
mation about the temporal history of the water produc- 
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lion from a single observation so that even the PVOUVS 
data gap (January 7 to January 30) is accessible to analy- 
sis. Moreover, the examination of the early March data 
can be undertaken on substantially smaller spatial scales 
(~2 x I0 6 km) than in the analysis of Stewart (1987). 

In this paper, the MCPTM corrected for optical depth 
effects is used to validate the ability of the model to accu- 
rately describe the distribution of hydrogen in the coma 
of Halley, to extract an improved set of H 2 O production 
rates for the PVOUVS observation period, and to obtain 
some information of the level of cometary activity during 
the period January 10 to January 30, where only OH 
radio measurements cover the entire interval. The obser- 
vational data are described in Section 2, followed by the 
modeling analysis in Section 3. In Section 4, model 
results are discussed, and conclusions from the analysis 
are presented. 

2. OBSERVATIONS AND DATA 

A comprehensive description of the observations and 
data is given in Stewart (1987). Here we shall summarize 
the most relevant points of that description. Comet Hal- 
ley was systematically observed by the PVOUVS at vari- 
ous wavelengths including Lyman-a from December 28, 

1985, to January 6, 1986, and from January 31 to March 7, 

1986. The PVOUVS instrument operated on a cycle syn- 
chronized with the spacecraft spin period of ~13 sec. 
During each cycle the instrument collected data through 
its 1.4° x 0.14° field of view as it scanned across the sky 
along a swath which varied from 20° to 80° in length. The 
data collected ideally represented the brightness of solar 
Lyman-a reasonantly scattered from atomic hydrogen in 
a swath of Comet Halley’s coma 1.4° wide and 20° to 80° 
long. The resolution along the swath was determined by 
the integration time which varied from 8 to 32 msec or 
0.19 10 to 0.75°. The data were collected for a period of 
1 hr and the resulting 230-270 swaths were averaged to- 
gether in order to improve signal-to-noise. A copy of this 
averaged data set, together with other ancillary informa- 
tion such as background, which consists of Lyman-a 
scattered by interplanetary hydrogen, was provided by 
A. I. F. Stewart. 

3. MODELING ANALYSIS 

The general approach adopted here is to simulate the 
spatial distribution of H atoms and Lyman-a resonance 
scattering in the coma of Halley and then compare this 
with the data described above. The simulation was ac- 
complished with the MCPTM which is described in detail 
in the earlier papers by Combi and Smyth (1988a,b; here- 
after referred to as Paper I and Paper II). For the parent 
molecule H 2 0 and its hydrogen-bearing daughter prod- 


ucts H and OH, MCPTM employs a spherically symmet- 
ric fluid model in the inner collisional coma interfaced 
with a Monte Carlo particle treatment for the transitional 
and collisionless regions of flow in the outer coma. The 
fluid model provides the density, velocity, and tempera- 
ture of the outflowing gas in the inner coma, which, in 
turn, are used by the Monte Carlo section of MCPTM to 
compute the velocity and spatial distribution of H in the 
outer coma. The distribution of H from H 2 0 and OH in 
the coma is determined by the photochemistry, elastic 
collisions, radiation force on H due to scattering of solar 
Lyman-a photons, the time-dependent production of 
H 2 0, and loss of H due to charge exchange with solar 
wind protons, ionization by solar wind electrons, and 
photoionization. While all these processes are described 
in detail in Paper I, the updating of certain rates and the 
fact that we are interested in Comet Halley rather than 
Comet Kohoutek require some brief rediscussion of the 
photochemistry, model production rate, and g factor. 

The kinetics of the photochemistry have changed 
slightly. The branching ratios for the various reactions 
are the same as in Paper II; however, the velocity distri- 
bution of the H and OH from H 2 0 dissociation has been 
altered in accordance with Crovisier (1989) who took into 
account the conversion of some of the dissociation en- 
ergy into vibration and rotation. The result is roughly a 
10% decrease in velocity of the fragments relative to that 
reported in Paper II. This modified photochemistry has 
also been used in all published MCPTM related work 
since Paper II. The composition of Halley’s atmosphere 
is principally H 2 0 (-80%). Other H-bearing molecules 
H 2 CO and HCOOH have been inferred to be present at 
levels < 10%, leaving H 2 O to generate >90% of the H. As 
a result, H 2 O was assumed to be the only source of H. 

In order to run the MCPTM, a production rate for H 2 O 
is required. Since the H 2 0 production rate varies sub- 
stantially over the several weeks needed to form the H 
coma, it is necessary to assume some time-dependent 
rate profile in the model. The rate profile shown in Fig. 1 
reflects the long-term variation of H 2 O production and 
was derived from a compilation of H 2 O production rates 
based on various Halley observations between 2.5 AU 
preperihelion to 2.5 AU postperihelion given in Combi 
(1989). This range easily encompasses the time period of 
interest here between late December of 1985 and early 
March of 1986. With the exception of the model calcula- 
tions shown by the solid line in Fig. 5 and by the triangles 
in Fig. 6, which will be discussed later, the production 
rate in Fig. I is employed in all model comae calcula- 
tions. In principle one should iterate the production rate 
history by utilizing the new history derived from compar- 
ing the data to the model comae to recalculate the 
MCPTM comae. The newly recalculated comae should 
then be compared to the data to derive the next iteration 
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70 


FIG. I. Model production rate. H : 0 production rate vs day of the 
year that is used in the MCPTM as derived from a variety of observa- 
tions and given by Combi (1989). 


of the history. This process should be repeated until con- 
vergence. In this study the iteration was terminated after 
the first iteration, since as shown in Fig. 2, the model and 
the data are already in close agreement. 

The g factor used in the MCPTM to calculate the 
brightness of the coma and the radiation force on the H 
atoms was taken to be given by 

photons atom' 1 sec -1 , 

where r(/) is the heliocentric distance in AU of the comet 
at time /, F(t) is the Lyman-a solar flux in units of pho- 
tons cm" 2 sec" 1 at the longitude appropriate to Comet 
Halley at time t y and f[v) is the Lyman-a flux profile 
(Lemaire et al. 1978) expressed as a function of heliocen- 
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FIG. 2. Comparison of model and data. The brightness profiles in kilorayleighs are shown along the swath as a function of the distance from the 
brightest point. The observations (* symbols) and corresponding MCPTM (solid lines) with a plane-parallel correction for optical thickness are 
compared for four profiles. The various figures (top left, top right, bottom left, bottom right) represent Jan. 4 (day 4), Feb. 9 (day 40), Feb. 24 (day 
55), and March 6 (day 65), respectively. The dashed lines show the background values, already subtracted from the data, for the Lyman-a scattered 
by the interplanetary hydrogen as provided by Stewart. These background values are typical of PVOUVS observations acquired on other days. It 
appears that there is still some additional background not included in the interplanetary hydrogen values. The Sun is to the right. 
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trie velocity and normalized so that its value at line cen- 
ter, i.e., the dimensionless quantity f(v = 0). is unity. The 
normalization of the 1.73 x 10"' f[v) factor corresponds 
to a total flux of 3.472 x 10" photons cm~ : sec"'. The 
term F(/)/3.472 x 10" effectively renormalizes gU) to 
the value of the Lyman-a flux present at the time of inter- 
est. F(/) for the period between December 1985 and April 
1986 was provided by Stewart (1991, private communica- 
tion), who used the 1989 compilation of SME (Solar 
Mesosphere Explorer) Lyman-a flux measurements to 
compute the fluxes at Halley’s longitude. 

Recently, Combi and Feldman (1993) have argued that 
the SME Lyman-a brightnesses should be —16% larger. 
They noted that observations by the SUSIM (Solar Ultra- 
violet Spectral Irradiance Monitor) experiment, which 
was designed explicitly to obtain accurate (better than 
5%) absolute calibration and was flown in August 1985, 
indicated a solar Lyman-a flux that was 16% greater than 
the value on the same date in the current version of the 
SME data base. Furthermore, their comparisons of 
nearly simultaneous IUE observations of OH and H in 
Comet Halley, and similar but previously published 
results for Comet P/Giacobini-Zinner, are brought into 
nearly exact agreement when the larger solar Lyman-a 
fluxes are used. As a result, 1.16F(/) was actually em- 
ployed in the calculations. Strictly speaking, one should 
employ the line center flux rather than the total Lyman-a 
flux in the formula for g. However, in the absence of 
measurements of the line center flux, using the total Ly- 
man-a flux is probably the best that can be done. More- 
over, the fact that the H atoms scatter Lyman-a photons 
anisotropically must also be taken into account. In the 
final calculation for the brightness of the H coma, the 
brightness output from the MCPTM was modified by 
multiplying by the following expression for the phase 
function of H atoms scattering Lyman-a photons, PU) = 
0.942 + 0.25 cos 2 (<K/)), where 6(t) is the Sun-comet- 
PVOUVS angle at the time, t, of observation. 

The particle trajectories for this model were integrated 
for a period of time equal to 4 x 10 6 sec, spanning a time 
interval from the time of observation to 4 x 10 6 sec be- 
fore. This is ample time for the hydrogen coma to fully 
develop. The number of particles employed in the simula- 
tion of a given swath was 20,000, which proved to be 
sufficient to provide favorable statistics where the data 
statistics were also favorable. The MCPTM grid was in 
each case a 6 x 10 5 km by 6 x 10 5 km square which for 
almost every measurement is smaller than the resolution 
along the swath. 

MCPTM generates a two-dimensional brightness dis- 
tribution of the solar Lyman-a flux scattered by the H 
coma in the skyplane of PVOUVS. As described above, 
the PVOUVS measures, however, a swath of brightness 
generated by a 1.4° x 0.14° slit scanning and integrating 


across the coma. Thus, the MCPTM output was con- 
volved with the PVOUVS slit function applied along the 
trajectory of the data swath to produce an MCPTM simu- 
lation of the PVOUVS scan. The model brightness along 
the swath was then normalized to the data by overlaying 
and shifting vertically until some optimal level of overlap 
occurred. The H 2 O production rate implied by the data 
was determined by the ratio of the normalized model 
profile to the unshifted model profile. Selected examples 
of this fitting procedure are shown in Fig. 2 and will be 
discussed later. This method of matching profiles appears 
to yield production rates which are reproducible to 
roughly 5%. Prior to matching the data, a background 
subtraction was performed. The Lyman-a background in- 
tensity which is primarily due to resonant scattering by 
the interplanetary hydrogen was provided by Stewart 
who used the Ajello (1990) model of Lyman-a scattered 
from this source. 

In practice, the profile fitting procedure was controlled 
principally by the data roughly between 1.5 x 10 6 and 3 x 
10 6 km. The brightest point and the points within 1.5 x 
10 6 km of that point were ignored in the overlaying proce- 
dure because of the uncertainty in the pointing direction 
of PVOUVS. (There is an uncertainty of about 0.25° in 
the position of the brightest point, i.e., the brightest point 
is not recorded at its true position in the coma.) For the 
other data points which are farther away, the conse- 
quences of pointing error are relatively less important 
since the brightness gradients diminish with distance 
from the nucleus. Points beyond 3 x 10 6 km were also 
considered in the profile fitting procedure, but were al- 
lowed less weight because of the increasing importance 
of noise and uncertainty in the amount of Lyman-a back- 
ground from the interplanetary hydrogen and because it 
was easier to decide the optimal overlap by including 
fewer points even when an excellent match was attaina- 
ble. 

Another potential problem with the analysis of the re- 
gion near the brightest point is that the hydrogen coma 
displays a notable optical thickness to Lyman-a within 
roughly 1.3 x 10* km of the brightest point. The proper 
treatment of this multiple scattering problem is beyond 
the scope of MCPTM and at least in principle is nontri- 
vial. As an attempt to contend with the multiple scatter- 
ing, an extra correction to the brightness generated by 
MCPTM was optionally included. The correction was 
based on the radiative transfer calculation of McElroy 
and Yung (1975). Their approach is valid for a semi-infi- 
nite plane parallel scattering atmosphere observed at nor- 
mal incidence which is not the same as the problem of 
scattering from the approximately spherical hydrogen 
coma of Halley observed at finite angles to the direction 
of insolation. Recently, however, a radiative transfer 
(RT) model, which was initially developed by Anderson 
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and Hord (1977), was applied to the hydrogen Lyman-a 
cometary coma by Combi and Feldman (1993). Although 
there are some significant assumptions in RT, such as a 
spherical coma and no bulk flow of H, near the nucleus 
where the optical depth is significant the RT cometary 
model is nevertheless quite accurate. The RT model was 
run for the conditions of February 24, 1986, and com- 
pared to the plane parallel corrected model used here. 
These results are discussed in Section 4. 

This profile fitting procedure was utilized on only a 
small fraction of the many Lyman-a swaths. The data 


TABLE I 

H 2 0 Production Rate for Comet Halley Determined from the 
Lyman-a Observations 


Date 


Heliocentric 

distance 

(AU) 

H : 0 production rate 
(KF 9 molecules sec" 1 ) 

1985 

Dec. 28 

1.06 

2.37 


29 

1.59 

2.41 


30 

1.56 

3.11 


31 

1.53 

2.86 

1986 

Jan. 1 

1.49 

2.78 


2 

1.46 

3.53 


3 

1.43 

3.71 


4 

1.39 

4.0 


5 

1.36 

4.03 


6 

1.32 

3.88 


Feb. 1 

L10 

11.3 


4 

0.59 

14.4 


8 

0.59 

14.9 


9 

0.59 

17.8 


10 

0.59 

17.9 


11 

0.59 

15.7 


12 

0.59 

15.7 


13 

0.59 

16.0 


14 

0.60 

18.0 


15 

0.60 

17.8 


16 

0.61 

17.1 


17 

0.61 

14.4 


18 

0.62 

12.5 


19 

0.63 

12.8 


20 

0.63 

14.2 


21 

0.64 

14.4 


22 

0.65 

14.3 


23 

0.66 

13.0 


24 

0.67 

14.0 


25 

0.68 

13.7 


26 

0.69 

13.6 


27 

0.70 

13.8 


28 

0.7! 

13.7 


March I 

0.73 

12.8 


2 

0.74 

11.2 


3 

0.75 

11.1 


4 

0.76 

9.49 


5 

0.78 

7.3 


6 

0.79 

8.01 


7 

0.81 

7.40 


that were either contaminated, for example, by light from 
Venus, or that passed too far from the nucleus, were 
discarded. The remaining data were divided into sets ac- 
cording to the day that they were measured. From these, 
only the brightest and longest swaths were compared to 
the appropriate model runs. Table I contains FLO pro- 
duction rates with one estimate per day for most of the 
days that PVOUVS observed the Lyman-a coma. 

4 . DISCUSSION AND CONCLUSION 

Figure 2 compares MCPTM model brightnesses (solid 
line) and the PVOUVS data brightnesses (asterisks) de- 
scribed above. Plotted is the log of the intensity against 
distance along the swath from the brightest point (the Sun 
is to the right). The plots (top left, top right, bottom left, 
bottom right) correspond to January 4, February 9, Feb- 
ruary 24, and March 6, 1986 (day of year 4, 40, 55, and 
65), respectively, and more or less span the time when 
PVOUVS was observing Halley’s hydrogen coma. Ex- 
amining the graphs it is clear that beyond about ± 10 mil- 
lion km from the brightest point noise dominates the data 
and a comparison is not possible. The wiggles in the 
model at similar distances are the statistical fluctuations 
arising from the Monte Carlo nature of the simulation. At 
smaller distances, the model matches the data very well 
except possibly for the brightest point and/or its immedi- 
ate neighborhood. 

The case in point is Fig. 2 top right, which corresponds 
to February 9. Here the second brightest point is too far 
to the inside of the model curve. Since this same pattern 
persists on neighboring days, it is unlikely to be an arti- 
fact. A possible explanation is that the data reflect a 
change in the production rate that is not contained in the 
production rate function (Fig. 1) adopted in the MCPTM. 
As mentioned in the previous section, the MCPTM effec- 
tively employs a secular production rate which repre- 
sents the long-term trend of H 2 0 emission from Halley. 
In reality the actual release of H 2 O is known to be vari- 
able on time scales of days or less. Since Halley was 
particularly close to PVOUVS on February 9 (0.39 AU; 
on Jan. 4, Halley was 1.39 AU from PVOUVS for com- 
parison), the slit was able to resolve smaller distances 
and thus was most sensitive to changes in the production 
rate on relatively short time scales. The other examples, 
particularly February 24, appear well matched by the 
model everywhere. 

Figure 3 displays the brightness along a swath on Feb- 
ruary 24 for three different models. February 24 was cho- 
sen because the coma was relatively optically thick on 
that date. The dashed line represents the RT model run 
for a production rate of 1.6 x 10 30 sec' 1 and a tempera- 
ture of 14,000 K for the hydrogen. Since the opacity de- 
pends on the velocity distribution of the hydrogen atoms, 
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FIG. 3. Effect of optical depth. The brightness profiles in kiloray- 
leighs are shown along the Feb. 24 swath as a function of the distance 
from the brightest point. The dotted line corresponds to the RT model 
while the lower solid line is the plane-parallel corrected MCPTM for 
optical depth, and the higher solid line is the MCPTM without any 
correction for optical depth. The Sun is to the right. 

a temperature is required by the radiative transfer calcu- 
lation to approximate the effect of the true velocity distri- 
bution. The temperature of 14,000 K was shown (Combi 
and Feldman 1992, 1993) to best represent the average H 
atom velocity distribution in the optically thick portion of 


the model. The sudden cutoff of the RT model beyond 
about 10 million km is simply due to reaching the end of 
the RT grid. The lower of the solid curves is the MCPTM 
model run for the same production rate and corrected 
using the results of McElroy and Yung (1975) with the 
same temperature. Clearly the corrected MCPTM and 
the RT model are quite similar in the region where the 
optical depth is significant. The shape of the profiles is 
the same, and the MCPTM is only slightly brighter (a few 
percent) than the RT model. Consequently, the corrected 
MCPTM is able to reasonably take into account the effect 
of multiple scattering by the hydrogen coma. It should 
not be assumed, however, that the corrected MCPTM 
replaces the RT model. The averaging effect resulting 
from the convolution with the PVOUVS slit is largely 
responsible for the agreement. The preconvolution two- 
dimensional model brightnesses are substantially less 
similar so that significant discrepancies would be ex- 
pected to arise in sufficiently high-resolution observa- 
tions of the optically thick region. Beyond the vicinity of 
the brightest point, the two models deviate due primarily 
to the fact that the RT model is a spherically symmetric 
model and the more realistic hydrogen distribution in the 
MCPTM is asymmetric because of radiation pressure ef- 
fects. The higher solid line in Fig. 3 is the uncorrected 
MCPTM model which assumes an optically thin atmo- 
sphere. The optically thin MCPTM is about 20% brighter 



FIG. 4. H 2 0 production rates for Comet Halley. The production rates in units of 10” molecules see -1 are shown as a function of the day of the 

year from Nov. 1, 1985, U> April 30, 1986, for a variety of measurements. The symbols refer to the following: (A) H 2 0 production rates determined 
from IUE data for OH emission by Combi et at. (1993), (□) H 2 0 production rates determined from the revised analysis (Smyth et at. 1993) of the 
6300-A emission data ferf «to m,c o*y8 cn originally published by Magee Sauer et at. (1990), (x) H 2 0 production rates determined from 18-cm OH 
radio observations by B^keWe-Morvan rr at. (1990), (O) HiO production rates determined from rocket observations of H Lyman-a emission by 
McCoy et at. (1992), (*) HjO production rates determined from 6300-A emission observations of atomic oxygen by Fink and DiSanti (1990), (•> 
H 2 0 production rates det crrrMnc< * m P apcr from MCPTM analysis of H Lyman-a observations in Table I. 
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than the MCPTM corrected for optical ddpth near the 
center. The correction will be less for March and early 
January since the activity of the comet is lower, leading 
to a thinner H coma. It will also be smaller since the 
comet PVOUVS distance is greater and thus the spatial 
resolution is less. Larger corrections —40% will occur for 
data taken near perihelion (Feb. 9) when the comet is 
highly active and the comet PVOUVS distance is 
smaller. In any case, the corrected MCPTM model ap- 
pears to be able to generate the proper brightnesses. 

In Fig. 4, the H 2 O production rates in Table I are com- 
pared with those rates obtained by other observers. Peri- 
helion for Comet Halley, which is on February 9 (day of 
year 40), is indicated by the dotted vertical line. The solid 
circles are for the H 2 0 production rates determined from 
the MCPTM analysis of the PVOUVS Lyman-a bright- 
ness measurements of Halley's hydrogen coma between 
late December 1985 and early March 1986. The other 
symbols represent H 2 0 production rates derived from 
other types of measurements for the period between 
early November 1985 and late April 1986. These do not 
exhaust all the measurements of Comet Halley's water 
production rates obtained in the relevant time interval, 
except possibly near perihelion. However, they are suffi- 
cient to illustrate how the new rates stand in reference to 
other determinations of the H 2 0 production rates. 

As is clear from Fig. 4, toward late December 1985 and 
early January 1986, our production rates blend in well 
with the IUE OH measurements (A symbols) by Combi et 
al. (1993), who also used the MCPTM approach. In early 
January 1986 there appears to be a substantial discrep- 
ancy between the revised H 2 0 analysis (□ symbols) of 
(Smyth et al . 1993) of the 6300-A emission data for 
atomic oxygen originally published by Magee-Sauer et al. 
(1990). However, the aperture employed by Magee-Sauer 
et al. was a 0.058° disk while the PVOUS aperture was a 
1.4° x 0.14° slit. Consequently, the sharp increase in the 
production rate obtained by (Smyth et al. 1993) would be 
expected to appear in the much larger PVOUVS aperture 
somewhat after the last PVOUVS observation on Janu- 
ary 6. Beyond January 6 there is a gap in the PVOUVS 
ftata lasting into January 31. Continuous acquisition of 
data resumed on February 1 (day of year 32). The 
MCPTM production rate rises above 1.1 x 10 30 sec -1 , 
until near March 1 (day of year 60), after which it falls 
and joins smoothly with the IUE-derived production 
rates of March 9 (day of year 68) and beyond. 

For the period between February 7 and March 7 (day of 
year 38 to 66), there are a few H 2 O production rate deter- 
minations from the rocket observations (O symbols) of 
Lyman-a of McCoy et al. (1992) and the O CD) ground- 
based observations (* symbols) of Fink and DiSanti 
(1990). Within their estimated errors of about 50%, these 
determinations are quite consistent with the MCPTM 


results. Finally, there is also an extensive set of 18-cm 
OH radio observations (x symbols) for this time interval 
recently converted to H;0 production rates by BockelCe- 
Morvan et al. (1990). These are in poor agreement with 
the MCPTM results during this time. The interpretation 
of OH radio data, however, is complicated by a number 
of difficulties with the most insidious being the quenching 
of OH excited states. This is most severe at high produc- 
tion rates. Indeed, Fig. 4 shows that the agreement be- 
tween MCPTM and the radio derived rates is only poor 
very near perihelion where the highest production rates 
occur. 

The production rates calculated from the OH radio 
data shown in Fig. 4 have a quenching correction in- 
cluded. Yet whenever the production rate is large, they 
still systematically tend to undercut the production rates 
derived from the other types of observations by a factor 
of about three. The calculation of quenching is based 
upon a scaling relation for the value of the quenching 
radius which is a function of the production rate of H 2 O 
(Schloerb 1988). It is assumed in the derivation of the 
scaling relation that the ratio of ions/neutrals, plasma 
velocity, and plasma temperature (velocity distribution) 
at the quenching radius do not depend on the production 
rate and more importantly, the quenching radius itself. 
This is quite probably incorrect, and an effective means of 
reducing OH 18-cm radio brightnesses to production 
rates must await a proper description of the cometary 
plasma for different production rates and heliocentric dis- 
tances. 



FIG. 5. Comparison of model to Feb. 8 data. The brightness profiles 
in kilorayleighs are shown along a Feb. 8 swath as a function of distance 
from brightest point. The observations are indicated by the « symbols. 
The MCPTM with a plane-parallel correction for optical thickness and 
for the production rate in Fig. I is indicated by the dotted line. The 
MCPTM with a production rate which is increased to 3.5 x 10“ mole- 
cules sec 1 between January 10 and January 30 is indicated by the solid 
line. The Sun is to the right. 
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If the factor of three correction is applied to the OH 
radio data around January 22, the result is that the sharp 
increase that the O('D) shows around January 10 to a 
value of almost 2 x 10 30 sec 1 on January 13 may actually 
continue to over 3 x I0 30 sec 1 by about January 22 and 
then fall to the PVOUVS rate on February l of 1.1 x I0 30 
sec -1 . Then Halley's activity would actually peak about 
21/2 weeks before perihelion, which contrasts with the 
usual expectation of maximum water emission when 
closest to the Sun. It must be emphasized, however, that 
any such extrapolation of the OH radio data should be 
regarded as speculative until proper means of interpret- 
ing the OH result are developed. There is, however, inde- 
pendent evidence in the PVOUVS measurements them- 
selves that the factor of 3 may be correct. 

Since the coma requires several weeks to develop due 
to its large spatial extent, the PVOUVS profiles must 
reflect the production rate history of the comet for a simi- 
lar period of time. Thus, the PVOUVS data nearest to the 
gap may be studied for evidence of the great increase in 
comet activity. Unfortunately, the data in this period 
seem to be rather noisy (see Fig. 2, top right) and also 
span a relatively small spatial scale. An exception at least 
with respect to the amount of noise is an observation 
taken on February 8 and shown in Fig. 5. The dotted line 
is the model with the usual underlying production rate 
(Fig. 1) which clearly underestimates the antisunward 
data beyond 10 7 km. The solid line has a production rate 
which is similar to that shown in Fig. 1 except that be- 
tween January 10 and January 30 a large production rate 
of 3.5 x 10 30 sec* 1 is assumed. Clearly, the match is 


significantly improved. Some caution must be exercised 
in interpreting Fig. 5. Although the background provided 
with the data was subtracted, it is apparent from Fig. 2 
that there still may be some additional background which 
is not taken into account. Therefore, it is likely that the 
data for February 8 in Fig. 5 are also somewhat high. As a 
result, the production rate of 3.5 x 10 30 sec 1 used in the 
model should be regarded as an upper bound. In order to 
obtain an excellent match between the MCPTM with the 
usual production rate and the data on the antisunward 
side, an extra background subtraction of 0.2 kR must be 
performed. Such a large additional correction (—50% of 
the initial background removed) results in an unsatisfac- 
tory fit on the sunward side, however, so that the back- 
ground cannot be adjusted to explain the increased tail- 
ward abundance in the coma. The presence of an 
elevated preperihelion production rate in the February 8 
profile is also supported by a preliminary reanalysis of 
the H Lyman-a image acquired by PVOUVS from Febru- 
ary 2 to February 6, 1986, which was analyzed earlier by 
Smyth et al. (1991) for the production rate in Fig. 1 . If an 
additional H 2 0 (triangular shaped) source with a time 
duration of ±6 days about a maximum production rate of 
3.6 x 10 30 sec -1 on January 23 is assumed, the newly 
computed Lyman-a contours provide a better fit to the 
observations, particularly in the tail portion of the coma, 
than the contours based on the H 2 0 production rate in 
Fig. 1. The presence of an elevated preperihelion produc- 
tion rate is therefore supported by the PVOUVS observa- 
tion. 

Figure 6 is a comparison between Stewart's (1987) H 2 0 



FIG. 6. Comparison of H 2 0 production rates for the MCPTM and Stewart. The solid dots represent the MCPTM Lyman-a results in Table I, 
and the asterisks are the earlier values of Stewart (1987). The MCPTM Lyman-a results for the production rate in Fig. 1 held constant after day of 
year 55 are shown by the A symbols. 
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production rates (asterisks) and those of the MCPTM 
(solid dots) which includes a correction for optical thick- 
ness. The MCPTM and Stewart results are fairly similar 
between about February 14 and March I (day of year 
45 and 60). In both cases, there is, on the average, a 
high production rate after perihelion which tends to dimin- 
ish somewhat after February 19 (day of year 50). 
Both models display oscillatory behavior. The MCPTM 
oscillations seem to be somewhat more subdued 
but seem to follow those of Stewart. The greatest 
differences lie at early and late times when the 
MCPTM is significantly higher and lower than Stewart, 
respectively. 

A comparison between the results in Fig. 6 is not 
straightforward because of the considerable differences 
not only in the models employed to convert data bright- 
nesses to production rates, but also in the manner of 
utilization of data for the derivation of production rates. 
The MCPTM attempts to treat the hydrogen realistically, 
accounting for both the dynamics and chemistry of the H 
atoms in a proper fashion as previously described. The 
Stewart model is a one-velocity (11.3 km sec -1 , the har- 
monic mean of the 8 and 18 km sec -1 principal velocity 
components of hydrogen from photodissociation, 
Stewart 1993, private communication) Haser model. The 
corrected MCPTM also accounts for optical depth ef- 
fects, the H scattering phase function, and the more re- 
cent SME (1989) solar Ly a fluxes plus a 16% further 
enhancement which results in a g factor nominally 10- 
20% higher than in Stewart. The treatment of the data is 
also dissimilar in important ways. Stewart (1987) aver- 
ages the data along the swath within ± 1 . 1 5° of the bright- 
est point. In addition, Stewart also combines the above 
averages from all the swaths which intercept the comet 
nucleus on a particular day. On the other hand in this 
study, while the individual swaths were considered, the 
brightest point and a neighborhood of about ±1.5 million 
km were ignored to minimize any potential errors intro- 
duced by the treatment of optical depth effects and the 
pointing uncertainties discussed earlier. 

For observations at early times (i.e., between January 
30 and February 14, before day of year 45), the lower 
rates inferred by Stewart (1987) are the result of examin- 
ing the region near the brightest point where the substan- 
tial optical depth shields some H from view. At later 
observational times, the higher rates deduced by Stewart 
may be a consequence of the swath averaging described 
above which corresponds to an area about the comet in- 
creasing to 5 million km x 8 million km by March 7 (day 
of year 66). This large area contains substantial contribu- 
tions of material produced by the nucleus at times before 
the day of interest when the production rate was rela- 
tively flat and hence may conceal the sharp decline sug- 
gested by the IUE observations (Fig. 4) which begins 


toward the end of the PVOUVS observation period. In 
this study, the production rate determination is con- 
trolled by points of a single swath at a few million kilome- 
ters from the intensity peak and may in March reflect 
more recently emitted material. Furthermore, the phase 
function is > I and also decreases the inferred production 
rate. It could, however, be reasonably pointed out that 
the sharp drop in the underlying production rate curve 
(Fig. I) may be responsible for the rapid decrease in the 
production rates beyond March 1 (day of year 60). 
In order to test this hypothesis, simulations were done 
with a production rate profile similar to that of Fig. 1 
except that beyond February 24 (day of year 55) the pro- 
duction rate was held constant. The triangles in Fig. 6 
are the result. Clearly a fall in H 2 0 production is still 
evident which indicates that the merging into the 
IUE observations is not an artifact of the choice 
of the underlying production rate curve used in the 
MCPTM. 

In conclusion, the simulation of PVOUVS measure- 
ments of Lyman-a scattered from the Halley coma is 
another example of the successful application of the 
MCPTM to cometary comae. The spatial brightness dis- 
tribution is very well reproduced as long as the underly- 
ing production rate is sufficiently representative of the 
actual activity of the nucleus. The production rates de- 
rived here are consistent with most of those from other 
types of measurements which overlap the PVOUVS data 
acquisition period. They also blend with the early March 
IUE results, in contrast to the earlier results of Stewart 
(1987). Moreover, if it is reasonable to suppose that an 
improved calculation of the quenching radius could alter 
the 18-cm OH radio results by a factor of ~2 to 3, then 
there is currently no reason to expect a fundamental dis- 
cordance with any other observations. 

The MCPTM rates, considering the substantial differ- 
ences in the modeling, corroborate Stewart (1987) rates 
between about February 10 and March 1, including the 7- 
day oscillation. In early February and in March the dif- 
ferences may be respectively reconciled by optical depth 
effects and differences in the spatial scale of the data- 
derived quantities which are compared to theory. Fi- 
nally, the spatially detailed analysis of February 8 
PVOUVS data suggests an extended period of high water 
production, possibly to almost twice the perihelion val- 
ues, for much of the time between January 10 and Janu- 
ary 30 when PVOUVS was not observing. This is consis- 
tent with the rapid rise seen in CH’D) data (Fig. 4) starting 
around January 10, the O('D) deduced water production 
rates of 2.90 x I0 30 and 2.68 x 10 30 sec -1 on January 16 
and January 17 respectively of Smyth el al. (1995), and 
production rates of over 3 x 10 30 sec -1 near January 22 if 
the OH radio derived H 2 0 production rates are low by a 
factor of —3. 
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